Statistics of the dissipated energy in driven single-electron transitions 
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We analyze the distribution of heat generated in driven single-electron transitions and discuss 
the related non-equilibrium work theorems. In the adiabatic limit, the heat distribution is shown 
to become Gaussian, with the heat noise that, in spite of thermal fluctuations, vanishes together 
with the average dissipated energy. We show that the transitions satisfy Jarzynski equality for 
arbitrary drive and calculate the probability of the negative heat values. We also derive a general 
condition on the heat distribution that generalizes the Bochkov-Kuzovlev equality and connects it 
to the Jarzynski equality. 



Landauer principle 0, [|[ , linking the erasure of infor- 
mation to the heat generated in a computation process, 
plays an important role both in understanding the foun- 
dations of thermodynamics (see, e.g., @ and references 
therein) and in practical attempts to realize thermody- 
namically reversible computation (see, e.g., Q). From 
physics point of view, a typical elementary step of a com- 
putation process is an externally-driven and controlled 
switching between two distinguishable states of a ther- 
modynamic system. From this perspective, one of the 
basic facts of statistical mechanics underlying the Lan- 
dauer principle is the statement that any system ma- 
nipulated adiabatically, at frequencies below its energy 
relaxation rate, remains in a state close to the instan- 
taneous thermal equilibrium. As a result, the energy is 
transferred into or out of the system reversibly, with the 
overall increase of the entropy "of the universe" that can 
be made arbitrarily small. Landauer principle addresses 
the average values, whereas the statistical fluctuations 
can lead to its violation in any given realization of the 
switching process |5|-[7|. These fluctuations are particu- 
larly important in small systems, making it interesting 
to understand not only the behavior of the averages, but 
also the fluctuations. The question of the distribution of 
the generated heat, or of the work done on the system, is 
addressed in fluctuation theorems that have been 

of interest in the non-equilibrium statistical mechanics 
at least for the last 30 years These theorems, how- 
ever, describe only some general characteristics of the 
heat /work distributions (discussed in more details be- 
low), while the distributions themselves are known only 
for some toy models (see, e.g., 12]). The aim of this 
work is to calculate the distribution of the generated 
heat for one of the most basic models of nanoscale in- 
formation processing: "single-electron tunneling" (SET) 
13]. There, the information is encoded in the position 
of individual electrons in a system of mesoscopic con- 
ductors, and is processed through electron tunneling be- 
tween the neighboring conductors. For the purpose of 
studying the heat distribution, this model combines sev- 
eral attractive features. It is understood to a high preci- 
sion, sufficient, e.g., for the development of metrological 
applications [lj, [l5|; it is considered as the basis of a 
practical scheme of reversible computing fl(il - [l8j . and is 



FIG. 1: Driven single-electron transitions. Left panel: elec- 
trons jump between two electrodes across a tunnel barrier, 
with an external potential controlling the energy bias AU of 
the transition. If an electron tunnels with the bias AU, this 
amount of energy is dissipated as heat. Right panel: a basic 
realization of the transitions is provided by the single-electron 
box, where the gate voltage V g creates the bias via capacitance 
C g . Electrons tunnel through the junction with capacitance 
C and create charge ne on the box island. 



sufficiently simple for explicit calculations. Our results 
show that, counter-intuitively, even for discrete electron 
transitions, the width of the generated heat distribution 
can be smaller than the energy fcgl 1 of thermal fluctua- 
tions, and vanishes for adiabatic switching together with 
the average dissipated energy. We also discuss the appli- 
cability to the SET switching of the " Jarzynski equality" 
(JE) [H, one of the better-known work theorems, and de- 
rive an equation that generalizes a physically appealing 
but more limite d eq uality derived earlier by Bochkov and 
Kuzovlev (BK) [ljf. 

Our discussion applies to driven single-electron transi- 
tions in general, although for the sake of clarity, we con- 
sider the specific example of a single-electron box (SEB) 
[Ifllifj, a small conducting island coupled through a tun- 
nel junction with conductance Gt and capacitance C to a 
large external electrode, and biased with respect to it by 
the gate voltage V g applied through the capacitance Cq, 
see Fig. [TJ The island can carry extra electric charge 
ne. In the simplest regime of relatively low tempera- 
tures T, /3" 1 <C £ c = e 2 /(2C s ), where /? = l/k B T and 
Cs = C + C g is the total capacitance of the box, and 
the gate voltage V g is restricted to the range < n g < 1, 
where n g = —C g V g /e, the box dynamics can be lim- 
ited to just two possible states, electron on the island, 
n = 1, and electron off the island, n = 0. We consider 



2 



the basic switching dynamics between these two charge 
states with the ramp of the gate voltage across its value 



Sn 



1/2 = at which these states are degenerate 



in energy. We first focus on the regime when the time 
dependence 8n g (t) is not very rapid on the large scale set 
by the energy Eq, so that the transition n = — > n = 1 
happens with certainty, while it can still be rapid on a 
smaller thermal scale /3 _1 . Such a transition is in general 
accompanied by dissipation of energy. The electrodes of 
the normal-metal SEB can be modeled as systems of non- 
interacting electrons. When an electron with energy e 
tunnels into or out of such a system, it carries an amount 
of heat e — /jL, where \x is the electron chemical potential 
of the electrode. The energy e is conserved in tunneling, 
and therefore the total amount of heat deposited in the 
two electrodes in electron transfer process is equal to the 
difference AU (t) between the chemical potentials of the 
electrodes at the time t of the transition. In the SEB, 
AU{t) — 2Ec8n g (t), and the total heat Q generated in 
one ramp of the gate voltage is 



Q = 2E C 



±Sn g ( n ) = 2E C J dtSn g (t)I(t) , (1) 



where the sum is over all back- and- forth (±) electron 
transitions between the states n = and n — 1 at the 
time instants n during one ramp, and I(t) = J2i ±£(it) 
is the associated particle current in the junction. The 
heat Q is an experimental observable, which can in prin- 
ciple be determined either by direct measurement, or by 
detecting the position of the gate n g (ri) at each tran- 
sition. To calculate the distribution of heat Q over 
an ensemble of transfer processes, we consider the sit- 
uation of small junction conductance, Gt *C e 2 /h, 
when the dynamics of tunneling can be described by 
Markovian master equation for the occupation proba- 
bilities p n of the states n, assuming first the tunneling 
rates T± = ±(G T /e 2 )A[//(l - e^ AU ) characteristic for 
the normal-metal junction [l3j . The master equation 
for the occupation probabilities p n , of the two charge 
states, n = 0,1 reduces then to p = Td — T^p for 
P = Pi - Pa, where = T + - T_ = GrAU/e 2 , and 
T s = T+ + T_ = (G T AU/e 2 ) coth(/3A£//2). It gives the 
probability pjk{t\, t^) for the system to be in the state k 
at time t 2 , if it was in the state j at time t\ < t 2 as 



(-l) k drT D {T)e-^ dT ' T ^"> 



(2) 



To obtain the quantitative results for distribution of the 
generated heat Q, we take the time dependence of the 
gate voltage in Eq. JT]) to be linear, 2Ec8n g {t) = -qt, on 
the relevant small scale set roughly by T. (It can still be 
non-linear on a larger scale set by Eq » Since the 

transitions, and hence the current I(t), are suppressed 
away from the resonance, t — ¥ ±oo, one can separate out 
the average heat (Q) and integrate by parts in Eq. ([1]), 



reducing the equation for heat fluctuations Q = Q — 
(Q) in this case to Q = r\ j dth(t), where h(t) is the 
random realization of the charge on the SEB relative to 
its average in the gate ramp, n(t) = I(t) — This 
means that the central moments of the heat fluctuations 
can be calculated through the correlation functions 
of the charge: 



(Q m )=V m / dt 1 ...dt m K<- m \t ll ...,t m ), 



(3) 



which in turn can be obtained directly from the solution 
@ of the master equation: 

K^ m \t u ...,t m ) =p i (*i)n i (ti)p i *(ii,i a )fifc(t2) 

■ ■■Pu(t m -i,t m )hi(t m ) , (4) 

where Pj(t) = paj (— oo, t) are the probabilities of the two 
charge states evolving from the state n = at t —> —oo, 
and fijit) = j — (n(t)) = j — pi(t) is the value of the 
charge fluctuation h(t) in the state j = 0, 1. Summation 
over all repeated indices is implied in Eq. (Q}. 

Equations © - (HJ) allow one, in principle, to find any 
central moment of the heat distribution. We calculate 
explicitly the heat noise <jq = (Q 2 ) 1 ^ 2 and the third 
moment Aq = (Q 3 ) 1 / 3 which coincides with the third 
cumulant of the distribution. For m — 2, one gets: 



KW(t 1 ,t 2 )=p a (t 1 )p 1 (t 1 )e-^ dTr * (T \ h 



< t 



2 ■ 



(5) 



Since the charge dynamics we are considering is essen- 
tially classical, the correlator at t\ > ty, is determined by 
the condition i^ 2 )(ii,t 2 ) = K {2 \t 2 ,ti). The correlator 
([5]) substituted in Eq. §3§ gives the heat noise oq. 

The first few moments of the generated heat are shown 
in Fig. [5] as functions of the dimensionless ramp rate 
v = rj(eP) 2 / Gt ■ Despite the thermal fluctuations of en- 
ergies of tunneling electrons, the heat noise <tq becomes 
smaller than fc^T and vanishes together with the average 
dissipated energy (Q) in the limit of adiabatic reversible 
evolution, v — > 0. This means that, similarly to continu- 
ous diffusive processes [13], for discrete SET transitions, 
the dissipation is suppressed in the adiabatic limit not 
only on average, but for the individual ramps as well. 
Quantitatively, the instantaneous equilibrium probabil- 
ity p = To /Ts obtained from Eq. ([2]) , combined with 
Eq. (J3J and (J5J) gives the heat noise for v < 1: 



o 2 Q = cv{k B T) 2 = 2k B T{Q), 



(6) 



with c = (8/ir 2 )J2n=o( 2n + - °- 85 - The sec- 

ond equality in (J6j) follows from comparison to the av- 
erage dissipated energy (Q) = r\ J dtt(h(t)) obtained 
using Eq. @ to find the first-order adiabatic correc- 
tion to (n(t)) = pi(t). This equality resembles the 
standard fluctuation-dissipation theorem relating (in this 
case, quasi-) equilibrium fluctuations of the dissipated 
energy and the average linear response to a slow change 
inducing this dissipation. Numerical evaluation of both 
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Eq. ([2]) gives after some algebra for m = 3: 



FIG. 2: The average and the higher moments of the distri- 
bution of the dissipated heat in single-electron transitions as 
functions of the ramp rate of the gate voltage driving the 
transitions. For comparison, open dots show the results (in- 
cluding the fourth central moment - dots without a line) of 
direct Monte Carlo simulations of the transitions. The simu- 
lations also give the overall shape of the distribution (upper 
inset). Increasing ramp rate, v — 0.1, 1, 10, leads to broaden- 
ing of the distribution. Probability of extracting energy from 
thermal fluctuations in the transitions is shown in the lower 
inset. Solid curve is the result of the Monte Carlo simulations, 
with the two asymptotics (dashed curves) given by Eq. (flQl) . 



quantities for arbitrary v (as in Fig. ^ shows that this 
relation holds very accurately up to v ~ 5. Used together 
with Eq. (TIB")) below, it also means that the distribution 
function p(Q) of the dissipated energy is Gaussian, 



p(Q) = G8/47r<Q)) 1 / 2 e-^W-<«» a / 4 W>. 



(7) 



This conclusion agrees with the behavior of Xq shown 
in Fig. [3J which, as described below, satisfies Xq oc v <C 
<jq oc i/v for v <C 1. In the adiabatic limit, Q is produced 
by many back-and-forth transitions and should indeed 
become Gaussian by the central limit theorem. 

In the opposite limit of the rapid gate voltage ramp, 
v ^> 1, the heat distribution can be found in the "single- 
jump" approximation assuming that in this case, the sys- 
tem has time for only one electron transition from the 
initial (n — 0) to the final (n = 1) charge state. Then, 



p(x = PQ) = 



xjv 
l-e~ 



3xp{-~ [ — 



-}• (8) 



One can check that this distribution reproduces well the 
large v behavior of the moments plotted in Fig. [5J The 
leading part of these asymptotes can be obtained analyt- 
ically by completely neglecting the 1/v parts of the dis- 
tribution, reducing it to p(Q) = (Q/a 2 )e~® l 2a , where 
a = \fyj '/? = {-qe 2 /Gt) 1 ^ 2 , and giving, e.g., for the third 
cumulant Xq = [(tt - 3)y^j2} 1/3 a. 

For arbitrary is, the third cumulant Xq can be calcu- 
lated from Eqs. Q and (fj|. Equation Q combined with 



Po(ti)pi(h)[pi(h 



Po(*2)] 
h<t 2 < t 3 . 



(9) 



In all other time intervals, the correlator is defined, 
similarly to K^ 2 \ by the condition that it is symmetric 
with respect to permutations of the time arguments. The 
integrals of Eq. © as required in Eq. © can be evaluated 
numerically and give Xq shown in Fig. [5] For v ^> 1, this 
can be done analytically, with the result agreeing with 
that of the large v approximation given above. In the 
adiabatic limit v -C 1, one can use the quasiequilibrium 
expression for the probabilities, po,i = (1 T ^d/^s)/2, 
and get Xq = c'vk B T, with d ~ 0.22. 

One of the interesting features of the statistics of 
the dissipated energy in the SET transitions consid- 
ered in this work is the existence of the Q < region, 
i.e. finite probability P of an SET transition taking 
place with extraction of energy from thermal fluctua- 
tions rather than with energy dissipation. This proba- 
bility could be the basis of implementation of an SET 
version of the "Maxwell's demon" (22j: One can envi- 
sion extracting electrical energy from thermal fluctua- 
tions in multi-junction SET circuits, with the help of a 
measurement/feedback loop, similar to the "information- 
to-energy conversion" demonstrated recently by applying 
rotating electric field to a dimeric particle [7| . The mag- 
nitude of this probability can be calculated for slow and 
rapid gate voltage ramps from Eqs. ([7]) and (|5J|: 



P 



_ j (l-erfKciz/S) 1 /*])^, f« 1, 
ir 2 /{6v), 



(10) 



In the adiabatic regime, the probability can be large, 
P — >• 1/2, at the cost of the typical value Q e of the ex- 
tracted energy being small on the scale of thermal energy, 
Q e ~ y^fc^T. In the opposite limit, v 3> 1, P is small, 
but Q e ~ ksT. For arbitrary u, P can be calculated 
by direct Monte Carlo simulation of the transitions, see 
Fig. [2] (lower inset), which shows good agreement with 
the asymptotic behavior (|10[) in the two limits. 

The SET transitions represent a very convenient sys- 
tem for testing the general work theorems @-[lll HI] of 
non-equilibrium statistical mechanics. To derive them, 
we adopt a more general set of assumptions than used 
above. We do not assume any specific energy dependence 
of the transition rates T± besides the detailed balance, 
r+( r = e /9AC/ ^ r ^r_(ri), and in principle, allow the tran- 
sitions to transfer more than one electron across the bar- 
rier. The following discussion is thus applicable, with 
small adjustments, also to normal metal/superconductor 
(NIS) junctions, which in the box geometry, in addition 
to SET transitions, demonstrate individual Andreev pro- 
cesses [24j transferring pairs of electrons. We also allow 
the gate voltage to have arbitrary time dependence, and 
the evolution of the SEB to start and end not in a defi- 
nite charge state. Instead of transition probabilities ©, 
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it is convenient now to represent Markovian evolution of 
the SEB in terms of all possible multi-transition paths 
between initial — T/2 and final +T/2 times, so that the 
average of any function f(Q) of heat is 

oo 1 

(/(<?)> =EEp»</(3W (ii) 

Here, as before, we restrict the dynamics to the two 
charge states, p n is the initial probability of the state 
n, and the averages (f(Q)) n ,j for trajectories starting at 
n and making j back-and-forth transitions are: 

/T/2 r-Tj ST2 

drj / dr 3 -_i... / drx 
-T/2 J-T/2 J -T/2 

Pn(n,T 2 , ...,Tj)f(Q(Ti, ...,Tj)). (12) 

The probability density p n can be expressed in terms of 
the probability f?±(Vi)7i+i) = e - -^ r±(t)dt ^.j^ 
system does not make a ± transition in the time interval 
from Tj to Tj-i-i. For odd values of j, 

P n (n, ...,Tj) = p ± (-r/2,Tl)/5 T (Tl,T 2 ) ■ ... • P±{Tj-l,Tj) 
Vt( T 3'T/2), p±{Ti,T i+1 ) = ^±(^,^+1)^(^+1), 

where the two sets of indices correspond to the two states 
n. For even j, the expression is constructed similarly. 

The detailed balance for tunneling rates implies im- 
mediately that the expression of the average of f(Q) = 
e -m _ jjj =i g^Atr(n) ] taken according to Eqs. ((TTJ> and 
(1121) , coincides with the transition probability w between 
the charge states in the time-reversed evolution [25| with 
the same voltage ramp, i.e. 

(e-"«> = £> (13) 
U 

where i, f — n(=pT/2) are the initial and final charge 
states of the system. This equation gives a general 
version of the work equalities discussed previously 0- 
ITlj |. m the form convenient for the SET transitions. In 
particular, if the gate voltage ramp satisfies the condi- 
tion 6n g (—t) = —5n g (t), as, e.g., does the linear ramp 
discussed above, the tunneling rate has the symmetry 
r_(— t) = T + (t). In this case, transition probabilities 



are symmetric, w(f — ► i) — w(i —> /), and Eq. (p~3|) 
gives (e~P®) = 1, the relation used above to derive the 
distribution ([7]), and by setting /3Q = AS, it is seen to 
be identical to the central equation (9.1) of Ref. [ll[ of 
Bochkov and Kuzovlev. We see that it is satisfied here 
for arbitrary distribution of the initial probabilities pi of 
the charge states. If the ramp is not antisymmetric rela- 
tive to the degeneracy point at t = 0, the right-hand-side 
of Eq. (|13p can be different from 1. 

For equilibrium initial state, the probabilities pi can be 
expressed through the equilibrium probabilities pf de- 
fined for the gate voltage at the end of the evolution: 
Pi = pfeP( Ef ~ Ei ~ AF \ where AF is the difference of the 
free energy for the final and the initial gate voltage values 
in the ramp, and Eij are the energies of the correspond- 
ing charge states given by the box charging Hamiltonian 
H(n) = Ecn 2 — U g n, where U g = eC g V g /C-£. This means 
that the energy difference between the initial and final 
charge states can be written as 

Ef-Ei = AH = - / \u g -E c )dn{t)- J ndU g (t), (14) 

where the first integral coincides with the heat ([1]), 
while the second term defines the " thermodynamic" work 
E HI Wth = - J ndUg. While this work is differ- 
ent from the actual physical work (2|| [27j done on the 
box by the source of the gate voltage, which is given by 
J Ugdn, it describes the change of the energies of the 
charge states. The average (e~ l3Wth ) over the distribu- 
tion of the thermodynamic work in a ramp can expressed 
similarly to Eq. (TT3l) . Relating then the probabilities pi 
through pf in this expression with the help of Eq. (|14l) , 
one gets the Jarzynski equality 

( e HS(Wth-AF)) = ! (15) 

for an arbitrary gate voltage ramp. Overheating effects 
were not considered in this analysis. 

In summary, we obtained the distribution of heat gen- 
erated in driven single-electron transitions and its central 
moments, and analyzed the most common work theorems 
and their validity in this context. The authors would 
like to acknowledge useful discussions with T. Ala-Nissila, 
K.K. Likharev, M. Mottonen and O.-P. Saira. 
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